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The volume W of the accessible N-body phase space and its dependence 



on the total energy is directly calculated. The famous Boltzmann relation 
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Q\ , S = k * ln{W) defines microcanonical thermodynamics (MT). We study how 



phase transitions appear in MT. Here we first develop the thermodynamics of 
microcanonical phase transitions of first and second order in systems which 
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O , are thermodynamically stable in the sense of van Hove. We show how both 

O 

kinds of phase transitions can unambiguously be identified in relatively small 
isolated systems of ~ 100 atoms by the shape of the microcanonical caloric 
equation of state < T(E/N) > and not so well by the coexistence of two 
spatially clearly separated phases. I.e. within microcanonical thermodynamics 
one does not need to go to the thermodynamic limit in order to identify phase 
transitions. In contrast to ordinary (canonical) thermodynamics of the bulk 
microcanonical thermodynamics (MT) gives an insight into the coexistence 
region. Here the form of the specific heat c(E/N) connects transitions of 
first and second order in a natural way. The essential three parameters which 
identify the transition to be of first order, the transition temperature Tt r , the 
latent heat qi a t, and the interphase surface entropy As sur f can very well be 
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determined in relatively small systems like clusters by MT. It turns out to 
be essential whether the cluster is studied canonically at constant tempera- 
ture or microcanonically at constant energy. Especially the study of phase 
separations like solid and liquid or, as studied here, liquid and gas is very 
natural in the microcanonical ensemble, whereas phase separations become 
exponentially suppressed within the canonical description. The phase tran- 
sition towards fragmentation is introduced. The general features of MT as 
applied to the fragmentation of atomic clusters are discussed. The similari- 
ties and differences to the boiling of macrosystems are pointed out. 

I. INTRODUCTION 

This is the first paper in a series of papers treating the topological structure of the 
N-body phase space of atomic clusters by microcanonical thermodynamics {MT). It was 
one of the primary issues of cluster physics to understand macrosystems starting from the 
properties and interactions of its constituents, the single atoms or molecules. Especially 
one likes to understand how the most spectacular changes in the thermodynamical behavior 
of macro systems, phase transitions, develop with increasing number of atoms from atomic 
dimensions towards the bulk. However, small clusters have also their specific features. They 
can be charged, their shape (surface) is important, they can fragment, they can rotate, etc. 

Studying clusters as thermodynamical systems we consider the phase space which can 
be reached by all constituents of the cluster. The topology of this total accessible N-body 
phase space reflects (or implies?) the behavior of many interacting systems because the 
dynamics is often chaotic. Then the evolution of many replica of the same system under 
identical macroscopic initial conditions follows the structure of the underlying N-body phase 
space. It is ergodic. In nuclear fragmentation the ergodicity is presumably due to the strong 
and short ranged friction between moving nuclei in close proximity. Friction between atomic 
clusters is yet unknown but probably it exists there also. 



2 



First we have to discuss the concept of thermodynamics of small systems in general and 
especially of phase transitions. It is important to realize that isolated clusters must be 
treated microcanonically. Usually there is no external heat- or particle bath which defines 
temperature, pressure, or chemical potential^. Micro canonical thermodynamics is the proper 
theory for isolated small systems. 

Phase transitions are well defined as singularities in various canonical or grandcanon- 
ical expectation values as function of e.g. the temperature in the thermodynamic limit 
liniAr^oo \n/v=q- E.g. the melting transition of bulk sodium at normal pressure shows up 
as a peak in the specific heat in the neighborhood of T m = 371K on top of which there is 
a qi at 5(T — T m ) singularity c.f. fig.[j]. With the normal finite resolution the delta-function 
manifests itself by a jump in the bulk-specific heat by qi at , the specific latent heat, see fig.|l]. 
There are different ways to identify the phase transiton of melting: First one realises that 
the systems changes its stucture, it melts. Another way is to measure e.g. the specific heat 
c(T) or the caloric equation of state < T(E/N) > and realizing that at the melting the 
temperature does not rise if more energy is put into the system up to the point where the 
whole system becomes liquid. I.e. that the caloric equation of state < T(E/N) > has an 
anomaly (here a plateau). 

The basic questions are now: Can one define phase transitions in finite, eventually small 
systems like atomic clusters. Is it possible to define phase transitions of small isolated sys- 
tems microcanonically as function of the energy not as function of the temperature? From 
which size on do phase transitions exist in such small systems? Is a cluster of some 50 to 100 
atoms big enough? We will show that it is in fact much easier to identify phase transitions in 



It is possible to prepare atomic clusters initially in a canonical ensemble by embedding them in 
an inert carrier gas Q. However, most probe-reactions on these clusters are too fast to keep the 
cluster during this time at constant temperature. Moreover frequent collisions with the carrier gas 
will distort the cluster so that one would not measure the properties of the cluster alone. 
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the microcanonical ensemble and that this is possible in astonishingly small systems. Small 
clusters can undergo a phase transition and one can unambiguously distinguish continuous 
(second order) from discontinuous (first order) transitions by the form of the microcanonical 
caloric equation of state < T(E/N) >. We will see that the main signal of a transition in fi- 
nite systems is not the singularity of some canonical expectation value in the thermodynamic 
limit but a specific anomaly in the smooth microcanonical equation of state < T(E/N) > 
which can be well identified already in pretty small systems. Of course, the transition may 
change its character when the number of particles increases. However, in the case of the 
10-states Potts model as well in the case of the boiling of sodium or potassium the bulk 
values of the three quantities above are clearly approached by the corresponding transitions 
in small systems/clusters of some 100 atoms. 

A microcanonical ensemble has some peculiarities: It does not have a positive definite 
heat capacity. In fact at a phase transition of first order the specific heat c(e) becomes 
negative in general. Therefore the classical signal of a peak in the specific heat is not useful 
to characterize a phase transition in small systems. 

The second peculiarity of small bound systems like atomic clusters is a new structural 
phase transition which does not exist in infinite homogeneous systems: At higher excitation 
they often do not simply boil into a gas of monomers. On the contrary, they may fragment 
into few relatively large fragments. In this case the size of the droplets at the transition point 
can be of the order of the size of the system itself e.g. in the case of fission. The fluctuation 
of the fragment sizes is substantial. Then one cannot ignore the droplets compared to the 
monatomic vapor anymore. The phase transition is not determined alone by the equilibrium 
of the homogeneous liquid with the homogeneous monatomic gas controlled by the equality 
of the chemical potentials of liquid and gas as we are used to in conventional (grandcanonical) 
thermodynamics. Often the number of droplets is similar to the number of monomers or 
even larger. Such a process was found recently in the decay of hot nuclei and is called 
multifragmentation. 

It is one of the aims of this series of papers to study the conditions for fragmentation 
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of hot and/or highly charged atomic clusters. However, also rapidly rotating clusters could 
fragment or even fission symmetrically. Therefore it is important to develop microcanonical 
thermodynamics including angular momentum as a fundamentally conserved quantity or 
order parameter. 

As today one is not familiar to microcanonical thermodynamics we will discuss in the 
following sections some of its basic features and show how phase coexistence and phase sep- 
arations show up in MT. Canonical and grandcanonical ensembles have the great advantage 
that one can treat many simple examples like systems which can be transformed into ideal 
noninteracting gases analytical. The prize one has to pay is that both ensembles average 
in the limit of large particle numbers over spatial inhomogeneities like e.g. phase separa- 
tion and interphase surfaces. Only under simple geometrical conditions this defect can be 
cured by the Maxwell construction. In contrast to this MT allows for such inhomogenuous 
situtions in a natural way. 

II. MICROCANONICAL THERMODYNAMICS 

Statistical mechanics as foundation of thermodynamics goes back to the works of Boltz- 
mann and of Gibbs. Boltzmann's gravestone has the famous epigraph: 

[S = k * InW] 

which is the most concise formulation of the essence of thermodynamics. W is the volume 
of the N-body phase space of the system (here in units of ((27ih) 3N ) or the partition sum of 
all possible quantum states of the N-body system: 

Qiv^HE^, (i) 

v 

and S is the N-body entropy of the system. 

Microcanonical thermodynamics explores the topology of the N-body phase space and 
determines how Qn depends on the fundamental globally conserved quantities namely total 
energy E = N * e, angular momentum L, mass (number of atoms) N, charge Z, linear 
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momentum p, and last not least the available spatial volume V of the system. This definition 
is the basic starting point of any thermodynamics since Boltzmann, see e.g. [||]. In his 
"Vorlesungen fiber Gastheorie, II. Teil" he calls the microcanonical ensemble the "Ergode" 
It is an entirely mechanistic definition saying that if we do not know anything more 
about a complicated interacting N-body system but the values of its globally conserved 
macroscopic parameters the probability to find it in a special phase space point (N-body 
quantum state) is uniform in the accessible phase space. 

The entropy is defined as the logarithm of Q see Boltzmann's epigraph 

S(E, V, N) = Ns(e = E/N) = ln(Q(E, V, N)), (2) 

we take Boltzmann's constant k = 1. The thermodynamic temperature T f ^ is defined by 

_ dS{E, V, N) _ ds(e) 

< J ~ dE " ~dT> [6) 

T tM = i. (4) 

By Laplace transform of Q(E, V, N) one steps from the "extensive" variables like E, V, N to 
the intensive ones like T, P, /i. E.g. the Gibbs grand partition function and the Gibbs grand 
potential are then 

z(p, p, fj) = JJJ o °° n(E, v, N ) e -^+pv-„N) dE dv dN ^ (5) 

G(P,P,fi) = -Tln[Z(P,P,fi)}. (6) 

I.e. Any sufficiently isolated part of the system has the same grand partition function 
G(T, P, /i), i.e. the system becomes translationally invariant or homogeneous in the bulk 
limit. In the same way one gets the canonical partition function and the Gibbs free energy 
as 

Z{(3, P, N) = JJ™ Q{E, V, N)e- pN{£+Pv) dE dV, (7) 
GF(T, P, N) = -T ln[Z(j3, P, N)}. (8) 
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Phase transitions in atomic clusters were quite early investigated microcanonically by 
molecular dynamical simulation of e.g. the melting transition in Ar 13 clusters see the pi- 
oneering work by Jellinek, Beck and Berry ||. In these early calculations the transition 
energy could be fixed by the anomaly of the "caloric" curve Ekin/N vs. E to t/N. Neither the 
entropy, nor the proper thermodynamic temperature, nor the order of the transition was 
determined. In contrast to these early work we will here present a proper calculation of the 
volume W resp. Qn of the accessible N-body phase space and its dependence on the total 
energy. We thus follow literally the prescription of Boltzmann's epigraph. This way, we will 
be able to determine all decisive parameters of the phase transition in question: T tr , qi at , 
a sur f and As feoi ;, the gain of entropy per evaporated atom, for the clusters. 

III. DIFFERENCES BETWEEN MICROCANONICAL AND CANONICAL 

ENSEMBLE 

According to van Hove a system of N particles interacting via short range two-body 
attractive forces with hard cores is thermodynamically stable, the thermodynamic limit 
of N, V — > oo\n/v=q exists for such systems, intensive quantities like the specific energy 
have finite limiting values |§ see also ||. Then, the thermodynamics derived from the 
microcanonical partition sum Q(E, N, V) and the one derived from the canonical Z(/3, P, N) 
or the grand canonical partition function Z(/3,P,fi) (usually) coincide. Outside of phase 
transitions of first order, the relative fluctuations AE/E, or Aq/q vanish oc 1/yN. 

This is naturally quite different for finite systems. However, even for bulk systems at 
phase transitions (T = T tT ) of first order do the microcanonical and the (grand) canonical 
ensemble differ essentially. In the (grand) canonical ensemble the energy fluctuations Ae per 
particle remain finite even in the thermodynamic limit ((Ae) 2 \T tr oc qf at , the specific latent 
heat). Consequently, the difference between the microcanonical and the (grand) canonical 
ensemble persists at transitions of first order and we must expect both ensembles to describe 
different physical situations. In section [TV] and especially in the subsection [IV A| we will 
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illuminate this important point in more detail. 

Systems interacting via long range forces like unscreened Coulomb repulsion between 
charged parts of equal sign, or the centrifugal force when they are rotating, or under the 
long range gravity, don't have a thermodynamic limit and must be described by the micro- 
canonical ensemble. Such systems fragment macroscopically into, in general, several regions 
of high density — condensed matter — and also into, in general, several regions of low 
density — vapor or may be empty space. Even though in the first two examples, Coulomb 
repulsion or repulsive centrifugal force, the systems are unstable they may disintegrate ini- 
tially slowly enough to establish a transient statistical equilibrium. This is evidently the 
case in the fragmentation of hot nuclei and we expect this to be quite similar in the case 
of fragmentation of atomic clusters. Differently from conventional thermodynamics where 
systems must be in a homogeneous phase at fixed temperature everywhere, here the sys- 
tem is most likely inhomogeneous. The inhomogeneities and their fluctuations are more 
important in characterizing the state of the system than any mean values. In contrast to 
thermodynamics of the homogeneous bulk, in small systems or large systems with long range 
forces the entropy connected to different partitions of the system is an important part of the 
total entropy. Familiar formulas like the one-particle entropy 



are useless for calculating the total entropy. 

IV. FIRST AND SECOND ORDER PHASE TRANSITIONS IN SMALL SYSTEMS 

Macroscopic systems have a discontinuity in the specific heat c^^T) at first order phase 
transitions. Cb u ik{T) may have a finite peak at T « T tr . On top of this there is a spike 
qi at 5(T — T tr ). With finite resolution it shows up as a jump in c^fe(T) by the latent heat 
qi a t. A typical example is the specific heat of bulk sodium at the melting transition, fig.|l]. 
In contrast, a transition of second order is continuous at the transition temperature where 
Cbuik{T) has (in the example of the Ising model) a logarithmic singularity in T — T tr . 




(9) 
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Microcanonical thermodynamics gives new and deep insight into this. It will further allow 
to extend the concept of phase transitions to systems not treated before by thermodynamics 
like systems with long range forces or strongly rotating systems. We begin with the discussion 
of microcanonical phase transitions in standard model systems in which phase transitions of 
first and second order are well known. As an example we take the 2-dimensional 10-states 
Potts model for which the asymptotic thermodynamics is even known analytically 0. The 
q-states Potts model is a generalization of the Ising model by allowing q instead of only 2 
spin values at each lattice point. In two dimensions it has a second order phase transition 
for q < 4 and a first order transition for q > 4 from a spin-ordered (in the following often 
called colloquially "liquid") to spin-disordered phase (called "gas"). 

In ref. |J we determined the three basic parameters of phase transitions of first order 
within microcanonical thermodynamics, the transition temperature T tr , the specific latent 
heat qiat, and the specific interphase surface entropy As sur f for the 10-states Potts model, 
i.e. a systems with a nearest neighbor coupling. It was demonstrated that for surprisingly 
small systems the values of these three parameters are closer to their asymptotic values than 
in the canonical ensemble. This is so because most of the finite-size scaling is due to the 
large, but trivial, exponent in the Laplace transform, eq.(^), from the microcanonical to the 
canonical partition sum PJlOl . 



From our studies of the Potts model || we made the plausible conjecture, which will serve 
us as working hypothesis: The two types of phase transitions in thermodynamically stable 
finite systems are distinguished by the form of the microcanonical caloric equation of state 
T t hd,(e): A transition of first order has a backbending caloric equation of state T -1 = 13(e), 
see figure ^b. In contrast a phase transition of second order has only a horizontal saddle 
point in < T(E/N) >. For a system with nearest neighbor interactions the area between 
(3(e) and the "Maxwell" line (3 tr = l/T tr is twice the interphase surface entropy As sur f 
H and the coefficient of surface tension is a sur f = Ti r A^ 1//3 As sur /. Following M.Fisher a 
nonvanishing surface tension qualifies the transition to be of first order in other words 
our experience with the Potts model suggests that the backbending (or S-bend) [|I^-[l4| 



is a 
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necessary signature for a phase transition of first order in a finite thermodynamically stable 
system. It is interesting to notice that the suppression of the configurations with phase 
coexistence compared to those with pure phases in the canonical ensemble is used since long 



to calculate the surface tension in general thermodynamical stable systems cf.e.g. ||1 5||: 



P(pUre) = Z(T) (10) 

P(mix) = — — (11) 

GfF m jx GrFj )ltre -|- GF sur f (12) 

N 2/3 a surf (13) 

-Ptr (mix) = e -N*As surf = e ~N 2 ^a aurf /T tr ^ ^ 

/ P tr (phasel ) P tr (phase2 ) 

where P(pure) is the probability to find the system at T tr in one of the pure phases, whereas 
P(mix) is the probability to find it in configurations with two separated phases. This 
generalizes to other systems our conjecture that the backbending of T(e) is a neccessary 
condition for the surface tension and for a transition of first order. Of course, there may 
also be other reasons for a backbending of T(e). 

The left darkened area is the defect of entropy As sur f = J e £2 (5de that the system 'pays' for 
introducing interphase surfaces, which it finally gets back when the whole system is converted 
to the new phase at e — e 3 , right darkened area, and the interphase surface disappears [||. 
In the bulk the transition is discontinuous as function of T or (3 and "jumps" from the liquid 
< £i) to the gas branch (e > e 3 ) of the caloric curve, fig.0. As a function of the specific 
energy e the transition is however continuous. With rising e the system passes smoothly 
from the liquid phase over a mixed phase with large coexisting fluctuations of the two phases 
( "gas bubbles" and "liquid droplets" ) towards the pure gas phase when the specific energy is 
increased by the specific latent heat qi a t- Figure ^|a shows the corresponding specific entropy 
s(e) = Jq f3{e')de' . Notice that the entropy is a monotonously rising function with rising 
energy in the region of energies shown. In figure §a we subtracted a linear function a + be in 
order to visualize the transition. Otherwise, on the scale of fig.^a one would not be able to 
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distinguish s(e) from a straight line. The transition is characterized by the convex intruder 
in s(e) of depth As sur f. The effect of the interphase surface tension is a somewhat slower 
increase of the entropy with rising specific energy as the system prefers to create interfaces 
between regions of different phases, spin-aligned and disordered ones. Figure ^|c shows the 
specific heat capacity 

C{£) = = WWe (15) 

as a function of the specific energy e. (Here numerical fluctuations in (3(e) in figure |2|b have 
been smoothed). One can see within the coexistence region of ei < e < e% (shaded area in 
the microcanonical specific heat has two poles and becomes even negative in between. 

The convex intruder in the specific entropy s(e), fig.^a, is forbidden by van Hove's the- 
orem in the canonical ensemble for an infinite number of particles || see also e.g. j3|. At 
energies in the region of the intruder, an infinite system is unstable against spontaneous 
devision into two parts ||. Conventional thermodynamics of the bulk is blind in this en- 
ergy interval. The Laplace transform eq.([7|) has no stationary point on the branch with the 
positive slope of (3(e). It can only see the branches of c(e) in the regions e < E\ and e > £3. 
Thus the canonical specific heat Cf, u ik(T) will be positive and approach finite values at the 
transition of first order. At T = T tr Cb u ik(T) has an additional peak = qi at 5(T — T tr ). 

In the canonical ensemble the phase transition is a sharp transition (the sharp 5(T — 
T tr ) in the specific heat) only in the thermodynamic limit. However, as a function of the 
specific energy the microcanonical ensemble transforms smoothly from one phase over the 
coexistence region into the other phase, see above, even in the thermodynamic limit. The 
nonvanishing width in energy of the phase coexistence persists even in the thermodynamic 
limit. Therefore it seems to me the singularity of c(T) reflects more a — sometimes even 
unpleasant — mathematical feature of the canonical ensemble to suppress inhomogeneities 
like phase coexistence than it does illuminate the underlying physics of the transition. This 
is completely contained in the (smooth) microcanonical equation of state. In fact it is exactly 
here where important information becomes hidden in the canonical ensemble. 
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Even though our conlusion from the Potts model is very convincing esp. in view of 
our remarks to eq.([T3|), it is not a general proof that the S-bend of the microcanonical 
caloric curve is neccessary for a phase transition of first order. S. Berry and others give 
some sufficient conditions for the backbending of the caloric curve in several model systems 
|TB| , P^JTB| -P^[ . These give a lot of insight into the mechanism that produces the S-bend of 
< T(E/N) >. Berry uses Nose-molecular dynamics to work with intensive variables like the 
temperature whereas we emphasize the advantages of the microcanonical ensemble and the 
use of the conserved extensive quantities like the total energy per particle e which are better 
suited especially for the study of phase separation. 

If the specific latent heat qi at — > and the specific interphase surface entropy As sur f — > 
the caloric equation of state gets only a saddle point at the transition. Then < E(T) > /N 
as well as < T(E/N = e) > become single valued, the transition is continuous in the 
canonical as well as in the microcanonical ensemble. We have a phase transition of second 
order. The two poles of c(e) merge and c(e) or Cb u ik(T) has a singularity at the transition 
point e tr ,T tr . Both slopes of c(e) are fully accessible in the canonical treatment of the heat 
capacity. Consequently, from the microcanonical caloric equation of state T(e) it is possible 
to identify and distinguish both kinds of transitions. In microcanonical thermodynamics the 
close relation between the two is very natural, transparent and simple. 

It is further instructive that in finite realizations of the two-dimensional Potts model with 
q = 10 spin orientations at each lattice point it was not possible to see a clean separation 
into a compact region of ordered spins and a compact region of disordered spins at energies 
inside the coexistence region even for a lattice of 100 * 100 points. There were always several 
"gas bubbles" and "droplets" fluctuating over the lattice and prohibiting large interphase 
surfaces. A clean interphase surface, another signal of a first order phase transition in bulk 
systems, cannot be seen in such small systems. Nevertheless the caloric equation of state 
T(e) is already close to its asymptotic form. (This is also the reason why it is possible to 
study many-body phenomena like phase transitions in relatively small systems with some 
100 atoms like we do here.) We can conclude from this observation that the other classical 
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signal of a transition of first order, a clear separation of the two phases, is not a useful signal 
of a transition of first order in small systems. 

A. Phase separation of macroscopic systems in the canonical ensemble 

According to eq.|T^ inhomogeneous configurations like configurations with two separated 
phases become exponentially suppressed by a factor exp(—Aa sur fN 2 ^ 3 /T) in the canonical 
ensemble. Nevertheless such configurations are frequently treated by conventional ther- 
modynamics with the help of the Maxwell-construction. This is a method to construct 
the microcanonical, not the canonical ensemble at energies where both phases coexist: Ho- 
mogeneous configurations with Nvi particles in the lower phase ("liquid") are linked to 
homogeneous configurations with Nv g particles in the energetically higher ( "gas" ) phase by 
the condition of equal chemical potentials fii = \x g with /x = TdGF/dN along a well de- 
fined interphase surface. After taking the thermodynamic limit N — > oo these homogeneous 
configurations can be obtained from the homogeneous canonical distributions at the lower 
(higher) temperature in the limit T — > T^. (T — > T^.). The mean total energy can be split 
into the energy of the homogeneous liquid with specific energy vi£\ and the homogeneous gas 
part with specific energy v g e g and vanishing fluctuations. Consequently the fluctuation of 
the total energy vanishes as well and the combined ensemble is the microcanonical ensemble 
with the specific total energy e = u t ei + v g e g . 

This is in contrast to a situation where the total system is treated as a single canonical 
ensemble without a forced phase separation. There the fluctuation of the total energy per 
particle remains proportional to the specific latent heat. In the latter case configurations 
with phase separation disappear. 

Clearly this method works only if an external field like the gravitational field orders the 
liquid below a given surface and the gas above. In a free cluster such a simple interphase 
surface can usually not be fixed and the above method to circumvent the exponential sup- 
pression of the inhomogeneous configuration with separated phases (liquid and gas or solid 
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and liquid) by the canonical ensemble is not possible. 



V. PHASE TRANSITION TOWARDS CLUSTER FRAGMENTATION 

There is a lot of interest in cluster fragmentation. Most of the work done considered the 
dynamical evolution of fragmentation. Starting with the atom-atom interaction one studies 
the explicit time-evolution of some small clusters with the help of Molecular Dynamics (MD) 
|i9H25| . We know from other systems that the dynamical evolution of a complicated many- 
body system is very much guided by the accessible phase-space. Then a statistical treatment 
explains the main outcome of such reactions. Moreover, a thermodynamic approach of 
cluster fragmentation offers the immediate connection to the thermodynamics of the bulk 
and especially to the most interesting thermodynamic phenomena: phase transitions. With 
finite clusters one can study the development of structural transitions like melting with 
increasing number of atoms see for example the illuminating review article by S. Berry in 
or the work of 



Most of these studies have been performed by embedding the cluster in a constant tem- 
perature heat bath. This is done in MD by using the Nose dynamics. From what was 
discussed above in the example of the Potts model the use of the canonical ensemble to 
describe phase transitions in finite clusters is not only unnatural but it seems especially in- 
convenient for discussing phase coexistence. We saw that inhomogeneous configuration are 
suppressed oc exp(—a sur fN 2 ^ 3 /T) relative to configurations with homogeneous phases in con- 
stant temperature realizations. Presumably this exponential suppression is the main reason 
for the vanishing of AT C = T m — Tf , the difference of the melting and freezing temperatures 
for large systems (N — > oo) discussed by Berry |16|]. Here microcanonical thermodynamics 
offers a promising and quite natural alternative to investigate this problem. 

Fragmentation is a new phase transition of finite systems. As other phase transitions of 
first order it is intimitely linked to the abilitiy of finite systems to become inhomogeneous 
[5Pfl. For infinitely many atoms it may become the liquid-gas transition. This will be the topic 
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of the fourth paper in this series 0. Here we are not so much interested in the well known 
physics of macroscopic phase-transitions but more in phase transitions and thermodynamics 
of finite systems. 

The thermodynamics of isolated metal clusters has many similarities with that of hot 
nuclei see ref. ]2"7|] . A new aspect of atomic cluster fragmentation, however, is the fact that 



charge and mass degrees of freedom are nearly decoupled. This is in sharp contrast to nuclear 
fragmentation where neutrons and protons are strongly coupled by the symmetry energy. 
It leads to a predominantly symmetric nuclear fission whereas for example doubly charged 
alkali clusters fission highly asymmetrically until strong shell effects of the daughters favor 
a more symmetrical split. The conditions for symmetric versus asymmetric fission will be 



discussed in the next paper of this series p8fl . A further peculiarity of clusters plays an 



important role in fragmentation: Atomic clusters have a much higher internal density of 
states than nuclei. The internal entropy Si n t can easily become larger than 10 per atom 
whereas in nuclei it is rarely Si nt > 2 per nucleon. This difference is of course due to the 
different statistics, fermionic statistics in a nucleus and bosonic statistics for the phonons 
in the cluster. Moreover, the internal structure of larger clusters has significant transitions 
like melting at which the internal spectrum shows strong enharmonicities. The bulk specific 
heat c p (T) is often considerably larger near the melting transition than the Debye limit for 
harmonic vibrations see figjl]. As shown in refs. p9| , |30| this additional internal entropy leads 



to an enhancement of fragmentation and also may enlarge the evaporation times by up to 
three orders of magnitude. 

As the charge and mass degrees of freedom are nearly decoupled in alkali clusters, it 
is possible to drive the system through the phase transition not only by increasing the 
excitation energy e but also by increasing its charge Z. Multiply charged but cold alkali 



clusters may Coulomb-explode, e.g. pT|j32]| . This is the topic of the third paper in this series 
|33|j . A third way to explore the fragmentation transition would be to subject the cluster to 
large rotations and disrupt it by the centrifugal force. All three conserved quantities, the 
energy E, the charge Z, and the angular momentum L are possible nonfluctuating order 
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parameters of this transition. At present it is experimentally not possible to transfer a large 
known angular momentum to a cluster. In the analogous case of nuclear fragmentation 



the disruptive effect of large rotations has just started to be discussed see e.g. |Q. In a 
couple of papers Jellinek and collaborators explored very early the effect of rotation on the 
isomerization and the behavior of clusters at low temperatures P5|- |T7| . Their interest 
aims into a similar direction as ours. They are concerned with the important problem how 
the low energetic properties of a cluster becomes modified by the rotation. We here want to 
study the effect of the two-dimensional stress by the centrifugal force on the fragmentation. 
It is further important to realize that the longrange centrifugal force induces inhomogeneities 
in the system which again demand a statistical treatment by microcanonical thermodynamics 
which leads to interesting pecularities as a kinetic temperature that can considerably exceed 
the thermodynamic temperature |34| . 



We sample the accessible N-body phase space by Mcrocanonical Metroplis Monte Carlo, 
MMMC . Its basic idea is as follows: Since every quantum state of the N-body system has 
an equal probability in the microcanonical ensemble one has to sample each quantum state. 
There is no hierarchy of a few important ones representing the main behavior of the system. 
In view of the tremendously large number of participating states this is impossible to do. 
Therefore it is essential to group the states by some common quality. In most experiments 
there are several observables which remain undetermined. In a typical atomic cluster ex- 
periment one counts the number of charged fragments but usually does not specify their 
internal excited state. Also one will usually not determine the momentum of each neutral 
evaporated monomer. Consequently, the natural sorting principle is the size of the unob- 
served phase space. As explained in great detail in refs. p9| , |38| , |39| MMMC samples only 
the relevant number of degrees of freedom. It treats the internal degrees of freedom of the 
fragments using the known internal specific heat. Each set of explicit degrees of freedom 
represents a whole subset of quantum states of the system and is weighted by the size of the 
corresponding part of the phase space. There are only a few millions of different configura- 
tions which are important for a specific reaction. MMMC gives a systematic method to find 
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them. It is thus possible to determine the size of the accessible N-body phase space i.e. the 
total entropy which is rather difficult to get e.g. by molecular dynamics. 

Let us now consider the fragmentation of sodium clusters with increasing number of 
atoms and charges. Figure |3] shows the average mass of the three largest fragments versus 
the specific excitation energy for the fragmentation of Nafoo- The thick solid curve gives 
the thermodynamic temperature in Kelvin (right scale). Below e = 0.5 eV/atom a singly 
charged big "evaporation" residue with dominating mass and additionally a small charged as 
well as several small neutral fragments are produced. Above 0.5eV/atom we find two small 
singly charged and several uncharged fragments of similar size (< 10) (multifragmentation). 
The caloric equation of state < T(E/N) > shows a significant backbending, which is a 
clear signal for a phase transition of first order at e ~ 0.5eV/atom. The two phases are 
characterized by the presence or absence of a dominating large fragment. 

Doubling the size of the cluster gives an even more pronounced phase transition, fig.^. 
The transition is shifted to slightly higher excitation energies. Evidently, fragmentation is 
easier in a smaller cluster with a large surface to volume ratio. The behavior of the mass 
distribution with excitation energy is very similar to the smaller cluster. In fig.|5] we show 
the fragmentation of an even larger cluster, Na^o- The mass distribution is like in the two 
previous cases, the transition is shifted further up to ~ 0.7eV/atom. 

Rising the charge from Z = 2+ to Z = 5+ shifts the transition point again down to e = 
0.55 eV/atom, fig.|| Rising the charge has a similar effect on the fragmentation transition 
as dropping the size of the cluster. All calculations were performed in the conducting sphere 
approximation to include the polarization of the charge distribution of Na n , see paper III 
|33| of this series for a detailed explanation of the method. 

Figure |7| shows a calculation using only the monopole-monopole (point charge) Coulomb 
interactions without the higher multipoles by the induced image charges. The transition 
is shifted towards slightly higher energies. The polarization (mirror charges) induces a 
net attraction. In contrast to similar calculations in ref. |25[] we used here the experimental 



binding energies for Na^ with q = 0,1+ and 2 < n < 21 given by ref. [{40"l , f4l| , whereas in |29| 
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we had only experimental values for neutral and charged dimers and trimers. This change 
induces a stronger upwards shift of the transition than the difference between conducting 
and point-charge distribution. The simple reason is the predominant production of small 
fragments with masses between n — 1 and n = 10, so that the transition is especially 
sensitive to small changes in the binding energies of these small fragments. 



VI. CONCLUSION 

Microcanonical thermodynamics is the basics of all thermodynamics, but only recently 
we came into the position to explore it for nontrivial situations. In contrast to the more 
familiar canonical thermodynamics it is the only formulation which is valid for small isolated 
systems like free atomic clusters. We can now study how phase transitions evolve in systems 
with rising number of atoms from small towards the macroscopic ones. 

Fragmentation phase transition is a new but generic phase transition of first order in 
small systems. Like other more familiar first order transitions as melting or boiling they 
lead to strong spatial inhomogeneities. However, unlike to the macroscopic transitions as 
boiling under the presence of an external gravity field, the inhomogeneities in microcanonical 
small systems are strongly fluctuating and the partition entropy plays an important role. 
Besides of being the correct ensemble for isolated systems the microcanonical ensemble is 
more suited to describe configurations with eventual phase separation than the canonical 
ensemble in which such configuration are exponentially suppressed oc exp(—N 2 ^ 3 a sur f /T). 
A fact that is especially important if one wants to study the evolution with rising particle 
numbers of configurations with two coexistent phase as e.g. the solid clusters with a molten 
surface or for the evolution of the liquid - gas transition as will be discussed in |2| . 

The fragmentation transition was first described for the multifragmentation of hot nuclei 
^ , |i3|j38| , [27[] . In this series of papers we will discuss the fragmentation phase transition in 



clusters under various external conditions. As we have seen this transition is well defined in 
sodium clusters of some 100 atoms. With rising size and increasing volume to surface ratio 
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it shifts towards higher excitation per atom. Charging the cluster leads to a lowering of the 
critical excitation energy per atom. 

In contrast to nuclear multifragmentation, the fragmentation of alkali clusters is very 
asymmetric if not special magic shell closures in the fragments favor more symmetric split- 
ting. This topic will be discussed in the next paper of this series . 

All our results concern the topological structure of the N-atom phase space only. Whether 
a realistic dynamical system can explore all the topological details of the phase space depends 
on how ergodic the dynamics is. The strong friction between nuclei moving one against 
another at close distances certainly helps to equilibrate the various degrees of freedom in 
nuclear fragmentation. Whether there is a similar strong friction between moving atomic 
clusters is not known. However, this is to be expected since atomic clusters have normally 
a less compact surface than nuclei. 

The basic limitations and approximations which we had to use in our analysis are dis- 



cussed in p7| and will be listed at relevant places in the following papers. 



VII. ACKNOWLEDGMENTS 

D.H.E.G. is very gratefull to A. Hiiller, S. GroBmann, S. Berry and J. Jellinek for en- 
lightning discussions about the basic principles of thermodynamics. We thank the Sonder- 
forschungsbereich SFB 337 of the Deutsche Forschungsgemeinschaft for substantial support. 
We are also grateful to the Freie Universitat Berlin for supporting us with computer time. 



19 



FIGURES 

FIG. 1. Specific heat of bulk sodium at atmospheric pressure from 44,iq|. The dashed line 
represents the specific heat calculated within the Debye model. T m is the melting and T v is the 
boiling transition. Notice the overshooting of c(T) over the Debye limit of harmonic vibrations at 
the melting point. 

FIG. 2. a) Specific entropy s(e) = J £ Pmicro(£)de versus the specific energy e for the 2-dim. 
Potts model with q = 10 on a 100* 100 lattice. In order to visualize the anomaly of the entropy the 
linear function a + be (a = 0.2119, b = 1.4185) was subtracted. Because we use periodic boundary 
conditions one needs two cuts to separate the phases and the depth of the convex intruder is twice 
the surface-entropy. 

b) Microcanonical caloric equation of state, (3 m i cro {e) = 1/T(e) as directly calculated by MMMC 

c) Specific heat c(e) = — (3 2 /(d(3/ds). The canonical ensemble of the bulk jumps over the shaded 
region between the vertical lines at e\ and £3. This is the region of the coexistence of two phases 
one with ordered spins, the other with disordered spins. Here c(e) has two poles and becomes 
negative in-between. The canonical thermodynamics is blind to this region. Notice that the poles 
are inside e\ < e < £3, i.e the canonical specific heat remains finite and positive as it should. 

FIG. 3. Average masses of the three largest fragments and the thermodynamic temperature T 



as function of the excitation energy per atom for the fragmentation of Na 2 



100 



FIG. 4. Same as figure ||] but for Na: 



2+ 

200- 



FIG. 5. Same as figure I but for Na; 



2+ 

500- 



FIG. 6. Same as figure ||] but for Na^oo- 
FIG. 7. Same as figure || but for Na^oo and point charge Coulomb interaction. 
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